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For a fermion gas with equally spaced energy levels that is subjected to a magnetic field, the 
particle density is calculated. The derivation is based on the path integral approach for identical 
particles, in combination with the inversion techniques for the generating function of the static 
response functions. Explicit results are presented for the ground state density as a function of the 
j^^ • magnetic field with a number of particles ranging from 1 to 45. 
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I. INTRODUCTION 



2 ' The explosive growth of mesoscopic physics in particular enabled to obtain a tunable number of electrons confined 

in semiconductor quantum dots p ] . Experiments on vertical quantum dots indicated that the confinement potential of 
' - single quantum dots is well described by a parabolic potential [gj . Although the electronic eigenstates and eigenvalues 
y ' are essentially given by the Fock-Darwin states B , detailed experimental data Bfl reveal the importance of correlation 
^ . effects on the ground state properties of the electrons in quantum dots. In order to take these correlation effects into 
^ ' account, various approximate theoretical methods HH] have been used. E.g., the eigenstates and eigenvalues of an 
'~~^ , harmonic interaction model including the effects of a magnetic field have been studied with operator techniques pOl , 
^^ • and the harmonic interaction model has been used to explain specific features in the addition spectrum of a quantum 
<^ [ dot in a magnetic field [Q. 

fvq • The harmonic interaction model is one of the rare examples for which the thermodynamical properties are exactly 

^^ [ soluble, including the boson or fermion statistics, in the presence of a magnetic field |I2]. It can also function as a trial 
^^ model for the variational treatment of systems with more realistic interactions with the aid of the Jensen-Feynman 

ON inequality p3| . An example of this approach can be found in Ref . p4| , where the Jensen-Feynman variational approach 
^^ is used to describe the Bosc-Einstcin condensation in a gas of ^^Rb and in a gas of ^Li atoms. The spin statistics of 
the harmonic interaction model |13] can be treated within the same many-body path integral formalism. 

In the present paper, we study the density of harmonically interacting electrons in a parabolic quantum dot in a 

j^ magnetic field, taking into account the electron correlation effects analytically. This analysis is a natural extension 

ri ^ of our previous investigation of the thermodynamical properties of a confined system of spin-polarized fermions in 

I . the presence of a magnetic fieldjO] , using a method that combines the path-integral formalism fl^ , the method of 

'^ ' symmetrized density matrices ]13|, and inversion techniques for genera ting functions |17|-p^ . Instead of using the 

Ch stochastic approach Ig^l with the Ito condition on the magnetic field [|Tj for calculating the path integral for N 

O identical interacting oscillators in a magnetic field, we relied on a detailed investigation of the classical equations of 

motion. The quantum mechanical corrections to their classical action are exactly taken into account. 

The model system of N harmonically interacting oscillators in a magnetic field is described by the Lagrangian (units 
with h and m equal to unity are used throughout this paper) 
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where lOc is the cyclotron frequency, and where the potential energy V results from a harmonic confinement potential 
and a harmonic two-body interaction 
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The two-body potential might be either attractive or repulsive depending on the plus sign or the minus sign in the 
two-body interaction. 

As a first step in the calculation of the path integral, we studied the case of distinguishable particles |l2] . Because 
the magnetic field only affects the equations of motion in the xy plane perpendicular to the magnetic field, the 
Lagrangian naturally decouples into two contributions L = L^y + L^ , where L^y contains the magnetic field. This 
allows to calculate the propagator in the xy-pl&ne independently from the propagator in the z-direction. Through 
a transformation to the center of mass reference frame, one obtains a set of N three-dimensional oscillators in a 
magnetic field. The center of mass is described by an oscillator with frequency fi. The remaining A^ — 1 oscillators, 
associated with the internal degrees of freedom, have a frequency 



= ^/Il^±NuP. 
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Note that the case of a repulsive interaction imposes a stability constraint on the confinement potential: f2 has to be 
large enough to keep the repelling particles together. The magnetic field does not affect these frequencies in the z 
direction. But in the xy plane the equations of motion reveal two renormalized frequencies due to the magnetic field 
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where loj^ — uJc/2 denotes the Larmor frequency. The index cm. refers to the renormalized frequency in the center of 
mass coordinates. 

Once the propagator for distinguishable particles is known, it is projected on the antisymmetric representation to 
obtain the fermion propagator |12| . 

This paper is organized as follows. In Sec. II the one-particle correlation function is calculated for identical 
harmonically interacting oscillators in a homogeneous magnetic field. In Sec. Ill special attention is paid to the 
fermion ground state density in the xy- and a:z-plane. In section IV some concluding remarks are given. 

II. STATIC RESPONSE PROPERTIES OF THE MODEL SYSTEM 

In the path-integral approach to quantum mechanics the expectation value of a an expression A (r) is given by 

Jdri---JdrNKi{f,f3\f,0)A{f) 
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where f is the 3A^-dimensional vector containing the coordinates ri, . . . , rjv of all the A^ particles, and Kj (r,/3| r', 0) 
denotes the propagator in the Euclidean time P — 1/ (kBT) , with ks denoting the Boltzmann constant and T 
the temperature. The subscript / emphasizes that identical particles (fermions or bosons) are considered. For the 
probability density and its Fourier transform this gives 
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Substituting the expression for the propagator Kj (r,/3| r, 0), one obtains the following integral for the Fourier trans- 
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where Zi {f3, N) — J dfKj (f, /3|r, 0) is the partition function of N identical particles, P denotes a permutation of the 
particle indices, and ^ — +1 for bosons and ^ = —1 for fermions. The propagators Kq, K^}, ^aji,,scm ^^^ ^ujl,s 
are known in closed form, and explicitly calculated in Ref. iQ. In order to obtain tractable expressions for nq, the 
summation over all possible permutations will be rewritten as a sum over all possible cycles. 

A. One particle expectation values 

For the one point correlation function a factor e"*""' has to be taken into account in each permutation when 
applying the cyclic decomposition of the permutations. Indicating the number of cycles of length ^ by M^, the cyclic 
decomposition for riq becomes 



NX, VNY, 



NZj (/?, N)J J {2^f K^ (^ ^z, /?! VNZ, o) /C.,. ( VNX, VNY, f3 



NX, VNY, 



j: y. iM,K, (k, q) i^^^^ (/c, (k))--^ n I , .,.M, (^- (k))-^' , (8) 

, ...M^, ^ £'^£ -"^^'- V<^ ^ 



with 



/C, (k, q) = y dr,+i . . . / dr^S (r,+i - n) e^i""! [] i^^^ , (x,+i, y^+i, /3| x,, y,, 0) K^ (z,+i, /3| Zj, 0) e-*"^"-^/^, (9) 

and /Cf (k, q = 0) = /C^ (k) is the same function as Eq. (17) of Ref. |12] in the calculation of the partition function 
Zj{P,N). We point out that the positive integers i and M^ (with 1 < i < N) have to satisfy the constraint 
J2(^Mi = N. Taking into account the semigroup property of the propagators K^^,s {xj+ItHj+ItP I ^jtUji^) £^nd 
K^ (zj+i, P\ Zj, 0), one recognizes in )Ci (k, q) the partition function (over a time interval ip) of a driven harmonic 
oscillator in a magnetic field 

ICiik,q)= f f f K^^,{x,y,ep\x,y,0)K^{z,ep\z,0)e- ^o" '^^-)'^^-)''- dxdydz , (10) 

with the driving force 

k ^^^ 
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This partition function is known in closed form and given by 

^ak, q) = 2 (cosh 1(3 s - cosh IPlul ) ^^^ \ 4s (cosh 1(3 s - cosh 1(3ojl ))\2 sinh ^ ^ ' / ' 
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Substituting the force fq (r) from (nil) into the above expressions for 4>q^^qy and 0^^ yields 
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The remaining integrations over k and R in n^ are Gaussian and easy to perform, eventually leading to 
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The exponential factor in n^ accounts for the center of mass contribution, and it becomes unity for non-interacting 
particles {w = fl). The factor fiq is the expectation value of J^i e"' '^' in the subspace of the relative coordinate system 
with its corresponding partition function Z/ (/3, N). 

We now introduce the generating function Qi (/3, u, q) = X]?^=o i^^i [Pi ^) '^q) "^ for tti^ Fourier transform of the 
density, as was done before ||l3,|l3l hi the absence of a magnetic field. 
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where ^/ (/3, u) = X]Ar=o ^^ (i^' ^) "" ^^ ^^^^ generating function for the partition function Z/ (/3, iV). After straight- 
forward algebra one is left with 
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It is noted that in the limit q ^ the sum rule fiq^a = 1 is indeed satisfied. The density n (r) in real space then 
becomes 
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The sum rule J dvn (r) — 1 for the density is easily verified. In the next subsection the fermion ground state density 
will be examined, by inverting the defining series for the generating function Qi (/3, u, q). Subsequently results will be 
presented for the ground state density in the a;j/-plane and in the xz-plane. 



III. GROUND STATE DENSITY 

Because of alternating signs in the recurrence relations for the partition functions Z/ (/3, N) , the equation (|22 
is not appropriate for numerical purposes, in particular in the low temperature limit. This sign problem can be 
circumvented by an alternative inversion of the generating function Qi using contour integration: 
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Substituting ( p^ ) into the above expression yields 
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with the short hand notations 

nl = ql/4s, 4^ql/4s, nl = ql/4w, b = e'^^', 6i=e-«^+"-), 62 = e-^(^---). (27) 

By expanding the Fourier transform of the density riq in powers of 5, 5i and 62 one arrives at 
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with (a) = r (a + p) /r (a) the Pochhammer symbol. The function ^jv (^) — e ^^^^j (/3, we*^) /S/ (/3, m) has previ- 
ously [Ol been obtained. Using (ph, the density (still at arbitrary temperature) then becomes 
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The ground state density at T = is then readily obtained by taking the Umit /? — > oo: 
"•T=o (r) 



^ V ^^^ exp (-as (x2 + y2) _ ^^,2) (32) 

oo oo oo oo oo f ^\^mi+m2+ni+n2 / i}\in 

X 



E E E E E ^""^ ^"'^ 



n n n ^^ ^ F (mi + 1) T (m2 + 1) T (m + 1) T (na + 1) T (m + 1) 

?n— ?ni—0 ni—0 7712— 712— "^ "^ 

y, y^ y^ (mi + Wj + l)fc^ (m,2 + ^2 + l)fc^ (m + l)fc 

'i^oi'oi'o r(fci + i)r(fc2 + i)r(fc + i) 

""^0,^0 ho (2(mi+m2) + l)_2,,r(gi + l)(2(ni+n2) + l)_2,,r(g2 + l)(2m+l)_2^r(g+l)' ^^^^ 

with 

TV TV 
o-=T7 — \ -, and d^— — --. (34) 

The sum rule / drnT=o (r) = 1 can easily be verified. 

A. Density in the xy-plane 

The ground state surface density n{x, y) = nT=a {x, y, z) /nT=o(0, 0, z) in the xy plane is cylindrically symmetric, 
which means that n{x = 0,j/) contains all the information. Figure 1 shows Nn{x = 0,y)/w for 1 up to 45 electrons 
for the non- interacting case {w = f2) at zero temperature. The factor 1/w makes this expression dimensionless (in 
the units with h — m = I used throughout this paper). The factor TV is introduced for clarity in the figure, in order 
to avoid too strongly overlapping curves. The densities for A^ = 1, 3, 6, 10, 15, 21, 28, 36, 45 electrons, i.e. for 
closed shell configurations, are emphasized by dashed lines. Oscillations appear in the density profiles, indicating 
concentric orbitals of increased density around the center of the confinement potential. The oscillations become more 
pronounced if the number of particles increases. Previously p2| the effects of an attractive (w > il) and a repulsive 
(w < fl) two-particle interaction have been studied in the absence of a magnetic field. It was seen that a repulsive 
interaction induced an expansion whereas an attractive interaction induced a contraction of the gas of interacting 
fermions. These effects also appear for non-zero magnetic fields, but for brevity reasons they are not plotted. 

With increasing magnetic field the oscillations in the density gradually become less pronounced and the confinement 
is enhanced as lu^ increases. The density profile changes whenever the magnetic susceptibility (which is proportional 
to dEc/diUL with Eg = X]£<£ ^ the ground state energy) exhibits a discontinuity as a function of the magnetic 
field. For luj^ ^ w the energy spectrum behaves like a Landau spectrum. This is illustrated for the non-interacting 
case {w = fl) in fig. 2a and fig. 2b for TV = 10 electrons, and in fig. 3a and fig. 3b for TV = 28 electrons. The insets 
show the magnetic susceptibility as a function of lul ■ 



B. Density in tiie xz-plaxie 

Similarly as for the density in the xy-plane, the density profile in the xz-plane changes whenever there is a jump in 
the magnetic susceptibility. Figure 4 shows density contour plots for 4 fermions {w — fl) for various magnetic fields. 
The asymmetry of the x- and z-direction due to the magnetic field along the z-axis is clearly revealed in the plots. 
In fig. 5 the density contours for 10 fermions (w = SI) in the xz-plane are illustrated. 

IV. CONCLUSION AND DISCUSSION 

In this paper we have presented analytical results for the density of spin-polarized harmonically interacting fermion 
oscillators in a magnetic field, taking the fermion statistics of the particles into account. The approach presented 
here is valid for any number of electrons and for any temperature. We concentrated on the ground state density for 
a number of particles ranging from 1 up to 45, although also higher particle numbers can be treated. The density, 
in the xy- as well as in the xz-plane, shows a magnetic field dependency that is governed by the discontinuity in 
the magnetic susceptibility. Oscillations are present in the density and they are more pronounced as the number of 
particles is increased. These oscillations are smoothed out and finally disappear with increasing Larmor frequency 
ujl ■ Whenever the magnetic field causes a discontinuity in the magnetic susceptibility, the density profile undergoes a 
sudden change, thus providing a means for characterizing the parameters of the system. To the best of our knowledge, 
the path integral approach used in this paper is the only method so far that provides this detailed information on the 
density for an interacting fermion system. 
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Figure captions 

Fig. 1: Scaled surface density Nn{x ~0,y)w for CI = w in the plane perpendicular to the magnetic field, as a 
function of the scaled distance y/yo (with yo ~ 1/ ^/w^ units with h = m = 1 are used) from the centre of the 
parabolic confinement potential for A'^ = 1, . . . , 45 fermions in the limit lol/w — > 0. The densities corresponding 
to 1, 3, 6, 10, 15, 21, 28, 36, 45 fermions (i.e. for closed shell configurations in the absence of a magnetic field) 
are indicated by dashed lines and the corresponding particle number is explicitly indicated. 

Fig. 2: a) Scaled surface density Nn [x — 0, y) /w for U = w in the plane perpendicular to the magnetic field, as a 
function of the scaled distance y/yo (with j/o = 1/v^) from the centre of the parabolic confinement potential 
for iV = 10 fermions and for various Larmor frequencies lol/w ~ 0.2, 0.6, 1.2, 3.0. The inset shows the magnetic 
susceptibility for A'^ = 10 fermions as a function of the Larmor frequency; b) Same as in Fig. 2a, but for 
lol/w = 0, 0.4, 0.9. 

Fig. 3: a) Scaled surface density Nn [x — 0, y) /w ior fl — w as a. function of the scaled distance y/yo (with 
yo = l/\/w) from the centre of the parabolic confinement potential for N = 28 fermions and for various Larmor 
frequencies lol/w — 0, 1.0, 2.0, 3.0. The inset shows the magnetic susceptibility for A^ = 28 fermions as a 
function of wl; b) Same as in Fig. 3a, but for lol/w ~ 0.5, 1.5, 2.5, 3.5. 

Fig. 4: Contourplots of Nn (x, z) /w for iV = 4 fermions and for 51 = w as a function of the scaled coordinates x/xq 
and z/zo with Xo = zo = 1/ ^/w, measured from the center of the parabolic confinement potential, for different 
values of the magnetic field: (a) lol/w = 0.2, (b) lol/w = 1.0, (c) lol/w = 1.5, and (d) lol/w = 3.0. 

Fig. 5: Contourplots of Nn {x, z) /w for A^ = 10 fermions and for il = w as a function of the scaled coordinates x/xo 
and z/ Zo with xo — zo — 1/ \/w, measured from the center of the parabolic confinement potential for different 
values of the magnetic field: (a) lol/w — 0.3, (b) lol/w = 0.5, (c) lol/w ~ 0.7, (d) lol/w = 1.2; (e) lol/w — 2.0, 
(f) ujl/w = 2.5, (g) ul/w = 4.0, (h) ujl/w = 5.0. 
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